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Abstract 


The production cross sections of top-quark pairs in association with massive vector bosons 
have been measured using data from pp collisions at aA = 8 TeV. The dataset corres¬ 
ponds to an integrated luminosity of 20.3 fb”' collected by the ATLAS detector in 2012 
at the LHC. Final states with two, three or four leptons are considered. A fit to the data 
considering the ttW and ttZ processes simultaneously yields a significance of S.Ocr (4.2cr) 
over fhe background-only hypofhesis for ttW (tiZ) producfion. The measured cross sections 
are cTttw = fb and cr„-z = ft). The background-only hypothesis with neither 

ttW nor tTZ production is excluded at V.lcr. All measurements are consistent with next-to- 
leading-order calculations for the ttW and ttZ processes. 
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1 Introduction 

The top quark is the heaviest known elementary particle, and its large coupling to the Higgs boson sug¬ 
gests that it might be closely connected to electroweak (EW) symmetry breaking. Despite the fact that 
the top quark was discovered two decades ago [1, 2] some of its properties, in particular, its coupling 
to the Z boson, have never been directly measured. Several extensions of the Standard Model, such 
as technicolour [3-7] or other scenarios with a strongly coupled Higgs sector [8] modify the top quark 
couplings. 

With the centre-of-mass energy and integrated luminosity of the collected data samples at the Large 
Hadron Collider (LHC), the processes in which the electroweak Standard Model bosons (y, Z, W and H) 
are produced in association with top quarks become experimentally accessible. Measurements of the 
ttZ, tty and ttH processes provide a means of directly determining top quark couplings to bosons [9- 
11], while the ttW process is a Standard Model (SM) source of same-sign dilepton events, which are a 
signature of many models of physics beyond the SM. Example leading-order Leynman diagrams for ttW 
and ttZ production at the LHC are shown in Ligure 1 . Previous searches for ttW and ttZ production at the 
LHC have been carried out by the CMS collaboration at = 7 TeV and = 8 TeV [12, 13]. 

This paper presents measurements of the ttW and ttZ cross sections based on an analysis of 20.3 fb“^ 
of proton-proton {pp) collision data at = 8 TeV collected by the ATLAS detector. Depending on 
the decays of the top quarks, W and Z bosons, between zero and four prompt, isolated leptons^ may be 


* In this note, lepton is used to denote electron or muon, including those coming from leptonic tau decays. 
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(a) (b) 

Figure 1: Example leading-order Feynman diagrams for (a) ttW and (b) ttZ production. 


produced. Channels with two (both with same-sign and opposite-sign charge), three, and four leptons are 
considered in this analysis. The opposite-sign (OS) dilepton, trilepton and tetralepton channels are mostly 
sensitive to ttZ production, while the same-sign (SS) dilepton channel targets ttW production. Table 1 
lists the analysis channels and the targeted decay modes of the ttW and ttZ processes. Each channel is 
divided into multiple analysis regions in order to enhance the sensitivity to the signal. A simultaneous fit 
is performed to all signal regions and selected control regions in the four channels to extract cross sections 
for ttW and ttZ production. 


Process 

tt decay 

Boson decay 

Channel 

1 

t 

N 


{i-vb){qqb) 

rv 

OS dilepton 

no 

+1 

{b^vbWvb) 

qq 

OS dilepton 

no 


{i-vb){qqb) 


SS dilepton 

no 


{t-vb){evb) 

t-v 

Trilepton 

no 


{e^vb){evb) 

qq 

OS dilepton 

no 

ttZ 

{qqb){qqb) 

er 

OS dilepton 

yes 

{€-vb){qqb) 

tr 

Trilepton 

yes 


{e-vb){evb) 

rr 

Tetralepton 

yes 


Table 1: List of ttW and ttZ decay modes and analysis channels targeting them. The last column indicates whether 
a final state lepton pair is expected from a Z boson decay. 


2 The ATLAS detector 

The ATLAS detector [14] consists of four main subsystems: an inner tracking system, electromagnetic 
(EM) and hadronic calorimeters, and a muon spectrometer. The inner detector provides tracking in- 
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formation from pixel and silicon microstrip detectors in the pseudorapidity^ range \r]\ <2.5 and from a 
transition radiation tracker (TRT) covering I 77 I < 2.0, all immersed in a 2 T magnetic field provided by a 
superconducting solenoid. The EM sampling calorimeter uses lead and liquid argon (LAr) and is divided 
into barrel (|? 7 | < 1.475) and endcap (1.375 < \r]\ < 3.2) regions. Hadron calorimetry is provided by a 
steel/scintillator-tile calorimeter, segmented into three barrel structures in the range I 77 I < 1.7, and two 
copper/LAr hadronic endcap calorimeters that cover the region 1.5 < I 77 I < 3.2. The solid angle cover¬ 
age is completed with forward copper/LAr and tungsten/LAr calorimeter modules, optimised for EM and 
hadronic measurements respectively, and covering the region 3.1 < \r]\ < 4.9. The muon spectrometer 
measures the deflection of muon tracks in the range It/I < 2.7 using multiple layers of high-precision 
tracking chambers located in toroidal magnetic fields of approximately 0.5 T and 1 T in the central and 
endcap regions of ATLAS, respectively. The muon spectrometer is also instrumented with separate trigger 
chambers covering \q\ < 2.4. 


3 Simulated event samples 

Monte Carlo (MC) samples are used to optimise the event selection and the choice of signal regions, 
and to model all signal and certain background processes. In the following, the simulation of signal 
and background events is described in detail. Eor all MC samples, the top quark mass is taken to be 
nit = 172.5 GeV, and the Higgs boson mass is set to 125 GeV. 

The ttV{V = W, Z) process is simulated using the Madgraph5 leading-order (LO) generator [15] with 
up to one additional parton, using the CTEQ6L1 [16] parton distribution function (PDE) set. Pythia 
6.425 [17] with the AUET2B underlying-event set of tunable parameters (tune) [18] is used to simu¬ 
late showering and hadronisation. The ttV samples are normalised to the inclusive next-to-leading-order 
(NLO) cross-section predictions, using Madgraph5_aMC@NLO [19], including the off-shell tiZfy* con¬ 
tribution and interference. An invariant mass of at least 5 GeV is required for any opposite-sign, same- 
flavour pair of leptons appearing in the matrix element. The obtained cross sections are CTttw = 232 + 32 fb 
and cTtiz = 215 ± 30 fb, compatible with other NLO QCD calculations [20, 21]. The quoted uncertainties 
include renormalisation and factorisation scale and PDE uncertainties, including as variations. 

The Alpgen v2.14 [22] LO generator and the CTEQ6L1 PDE set are used to simulate W/Z production. 
Parton showers and hadronisation are modelled with Pythia 6.425. The WjZ samples are generated 
with up to five additional light partons, separately for fP/Z, WIZ+bb, WjZ+cc and Wc, and normal¬ 
ised to the respective inclusive next-to-next-to-leading-order (NNLO) theoretical cross sections [23]. To 
avoid double-counting of partonic configurations generated by both the matrix-element calculation and 
the parton-shower evolution, a parton-jet matching scheme (MLM matching) [24] is employed. The 
overlap between WjZ+QQ (Q = b,c) events generated from the matrix-element calculation and those 
generated from parton-shower evolution in the fP/Z-i-light-jet samples is avoided via an algorithm based 
on the distance in rj-ip space between the heavy quarks: if AR{Q, Q) > 0.4, the matrix-element prediction 
is used, otherwise the parton-shower prediction is used. 


^ ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector 
and the z-axis coinciding with the axis of the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and 
the (/-axis points upward. Cylindrical coordinates (r,<p) are used in the transverse plane, tp being the azimuthal angle around 
the beam pipe. The pseudorapidity is defined in terms of the polar angle 6 as 77 = - In tan(6/2), and the distance between two 
objects ’mrj — (f) space is measured in terms of AR h ■\j{l^ri)^ + 
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Diboson samples are generated using the Sherpa 1.4.1 [25] generator with the CTIO PDF set [26], with 
massive b- and c-quarks and with up to three additional partons in the LO matrix element. Samples are 
normalised to their NLO QCD theoretical cross sections [27]. Alternative models for the diboson back¬ 
ground are provided by the POWHEG-BOX 2.0 [28-30] generator, which implements the NLO matrix 
elements, interfaced with Pythia 6.425 or Pythia 8.1 [31]. 

Simulated tt and single-top-quark backgrounds corresponding to the t-channel, Wt and 5-channel produc¬ 
tion mechanisms are generated using the POWHEG-BOX generator, with the CTIO PDE set. All samples 
are interfaced with Pythia 6.425 with the CTEQ6L1 PDE set and the Perugia2011C [32] underlying event 
tune. Overlaps between the tt and Wt final states are removed through the diagram removal scheme [33]. 
The tt sample is normalised to the Top-i--i-2.0 [34] theoretical calculation performed at NNLO in QCD 
that includes resummation of next-to-next-to-leading logarithmic soft gluon terms [35-39]. The single¬ 
top-quark samples are normalised to the approximate NNLO theoretical cross sections [40-42] calculated 
using the MSTW2008 NNLO PDE set [43, 44]. 

The production of a single top quark in association with a Z boson through the t- and 5-channels, of 
the WtZ process, and of a top quark pair in association with a W boson pair {tt WW) are simulated 
with Madgraph5 LO and the CTEQ6L1 PDE set. Madgraph is interfaced with Pythia 6.425 using the 
AUET2B tune and the CTEQ6L1 PDE set. The relevant samples are normalised to the NLO theoretical 
predictions calculated with Madgraph5_aMC@NLO [19]. The production of three vector bosons that 
decay to three or four leptons is also simulated with Madgraph5 and Pythia 6.425. The LO cross section 
obtained from the generator is used to normalise the samples. The production of two W bosons with 
the same charge is modelled using the Sherpa generator, including diagrams of order and a|^Qfg. 
The LO cross section obtained from the generator is used to normalise the samples. The four-top-quark 
process (tttt) is simulated with Madgraph5 interfaced with Pythia 8. 

Associated ttH production is simulated using NLO matrix elements obtained from the HELAC-Oneloop 
package [45]. The POWHEG-BOX program served as an interface for shower MC programs. Samples 
were produced using the CTIONLO PDE set and showered with Pythia 8.1 with the CTEQ6L1 PDE and 
the AU2 underlying-event tune [46]. The ttH cross section and Higgs boson decay branching fractions are 
taken from the theoretical calculations collected in Ref. [47]. The process gg ^ H ^ 4b is modelled using 
the POWHEG-BOX interfaced with Pythia 8.1. WH and ZH production are modelled using Pythia 8.1. 
The samples are normalised to the NNLO QCD cross sections with NLO electroweak corrections [47]. 

All simulated samples produced with Pythia use Photos 2.15 [48] to simulate photon radiation and 
Tauola 1.20 [49] to simulate t decays. Events from minimum-bias interactions from the same bunch 
crossing as the hard-scattering process and in neighbouring bunch crossings, known as pile-up, are simu¬ 
lated with the Pythia 8.1 generator with the MSTW2008 LO PDE set and the AUET2 [50] tune. These are 
superimposed on the simulated hard-scatter events in a manner which reproduces the luminosity profile 
of the recorded data. 

All samples are processed through a simulation of the detector geometry and response [51] either using 
Geant4 [52], or Geant4 with a fast simulation of the calorimeter response [53]. All samples are pro¬ 
cessed by the same reconstruction software as the data. Simulated events are corrected so that the object 
identification, reconstruction and trigger efficiencies, energy scales and energy resolutions match those 
determined from data control samples. 
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4 Object reconstruction 


The final states of interest in this analysis contain electrons, muons, jets, t^-jets and missing transverse 
momentum. 

Electron candidates [54] are reconstructed from energy deposits (clusters) in the EM calorimeter that are 
associated with reconstructed tracks in the inner detector. The electrons are required to have |?7ciusterl < 
2.47, where /yduster is the pseudorapidity of the calorimeter energy deposit associated with the electron 
candidate. Candidates in the EM calorimeter barrel/endcap transition region 1.37 < |?7ciusterl < 1-52 are 
excluded. The electron identification relies on a likelihood-based selection [55]. 

To reduce the background from misidentified or non-prompt (labelled as “fake” throughout this paper) 
electrons, i.e. from decays of hadrons (including heavy flavour), electron candidates are required to be 
isolated. In the opposite-sign dilepton and tetralepton channels, in which such background is small, the 
electron isolation is defined using only tracking information. In the opposite-sign dilepton channel, the 
ratio of the sum of track transverse momenta in a cone of size AR = 0.3 around the electron 

track, excluding the electron track itself, to the transverse momentum (p!j) of the electron is required to 
be less than 0.12. In the tetralepton channel, the requirement is loosened to < 0.18. In the 

trilepton and same-sign dilepton channels, in which the background with fake leptons is more prominent, 
additional requirements are imposed on the electron isolation. Eor electrons with < 50 GeV, both the 
ratio of the additional calorimeter energy within a cone of AR = 0.2 around the electron (E.^^® ^) to the 
of the electron, and IPj are required to be less than 0.12. Eor pj > 50 GeV, both and 

p^<03 required to be less than 6 GeV. 

Muon candidates are reconstructed from track segments in the various layers of the muon spectrometer, 
and matched with tracks identified in the inner detector [56]. The final muon candidates are refitted 
using the complete track information from both detector systems, and are required to have \t]\ < 2.5. 
Additionally, muons are required to be separated by AR > 0.4 from any jet and to satisfy a pT-dependent 
track-based isolation requirement [57] that has good performance under high pile-up conditions. This 
requires that the scalar sum of the track transverse momenta in a cone of variable size AR = (10 GeV/ p^) 
around the muon (excluding the muon track itself) must be less than 0.05pf^. 

Eor both the electrons and muons, the track longitudinal impact parameter with respect to the primary 
vertex,^ ^o^ is required to be less than 2 mm. In the same-sign dilepton channel, in which backgrounds 
from fake leptons are dominant, it is also required to satisfy |zo sin 6\ < 0.4 mm, and the significance of 
the transverse impact parameter is required to satisfy \dQlo-{dQ)\ < 3, where cr{dQ) is the uncertainty on 
do- 

Jets are reconstructed with the anti-kf algorithm [58-60] with radius parameter R = 0.4 from calibrated 
topological clusters [14] built from energy deposits in the calorimeters. Prior to jet finding, a local cluster 
calibration scheme [61, 62] is applied to correct the topological cluster energies for the effects of non¬ 
compensating calorimeter response, dead material and out-of-cluster leakage. The jets are calibrated to 
restore the jet energy scale to that of jets reconstructed from stable simulated particles, using energy- and 
p-dependent calibration factors derived from simulations. Additional corrections to account for residual 
differences between simulation and data are applied [63]. After calibration, jets are required to have 
Pt > 25 GeV and |p| < 2.5. To avoid selecting jets from pile-up interactions, an additional requirement, 

^ A primary vertex candidate is defined as a vertex with at least five associated tracks, consistent with the beam collision 
region. If more than one such vertex is found, the vertex candidate with the largest sum of squared transverse momenta of its 
associated tracks is taken as the primary vertex. 
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referred to as the jet vertex fraction criterion (JVF), is imposed on jets with pT < 50 GeV and I 77 I < 2.4. 
It requires that at least 50% of the scalar sum of the transverse momenta of tracks with px > 1 GeV, 
associated with a jet, comes from tracks compatible with originating from the primary vertex. During 
jet reconstruction, no distinction is made between identified electrons and jet energy deposits. Therefore, 
if any of the jets lie AR < 0.2 from an electron, the closest jet is discarded in order to avoid double 
counting of electrons as jets. After this overlap removal, electrons and muons which lie AR < 0.4 from 
any remaining jet are removed. 

Jets containing fj-hadrons are tagged by an algorithm (MVl) [64] that uses multivariate techniques to 
compute weights by combining information from the impact parameters of displaced tracks as well as to¬ 
pological properties of secondary and tertiary decay vertices reconstructed within the jet. Larger weights 
indicate that a jet is more likely to contain fj-hadrons. The working point used for this measurement cor¬ 
responds to 70% efficiency fo fag a fj-quark jef, as defermined for fj-jefs wifh px > 20 GeV and \t]\ < 2.5 
in simulafed tt evenfs. The rejecfion factors for lighf-jefs and c-quark jefs are approximafely 130 and 
5, respectively. The efficiency of f^-fagging in simulafion is corrected fo fhaf in dafa using a tt based 
calibration [65]. 

The missing fransverse momenfum wifh magnifude is reconsfrucfed [ 66 ] as fhe negafive 

sum of fransverse momenfa of all elecfrons, muons, jefs and calibrated calorimefer energy clusters nol 
associafed wifh any of fhese objecfs. 


5 Event selection and background estimation 

The measuremenfs presented here are based on dafa collected by fhe ATLAS experimenf in pp collisions 
af = 8 TeV in 2012. The corresponding integrated luminosity is 20.3 fb“^ Only evenfs collected 
using a single-elecfron or single-muon frigger under sfable beam condifions, fhaf safisfy fhe sfandard dafa 
qualify criferia, are accepted. The frigger px thresholds are 24 or 60 GeV for electrons and 24 or 36 GeV 
for muons: the triggers with the lower px thresholds include isolation requirements on the candidate 
lepton, resulting in inefficiencies af high px fhaf are recovered by fhe friggers wifh higher px fhresholds. 
Evenfs are required fo have af leasf one reconsfrucfed primary vertex. In all selecfions considered, af 
leasf one lepfon wifh px > 25 GeV is required fo mafch (AR < 0.15) a lepfon wifh fhe same flavour, 
reconsfrucfed by fhe frigger algorifhm. 

Four channels are defined based on fhe number and charges of fhe reconsfrucfed leptons, which are sorted 
according fo fheir fransverse momenfum in decreasing order. For fhe opposife-sign dilepfon channel, 
fwo lepfons wifh opposife charge and px > 15 GeV are required. In fhe same-sign dilepfon channel, 
evenfs are required fo confain fwo same-sign lepfons wifh px > 25 GeV. In bofh dilepfon channels, 
evenfs confaining addifional lepfons wifh px > 15 GeV are rejected. For fhe frilepfon channel, evenfs are 
required fo confain fhree lepfons wifh px > 15 GeV. For fhe fefralepfon channel, exacfly four lepfons 
wifh Px > 7 GeV are required. Evenfs safisfying bofh fhe frilepfon and fefralepfon channel selecfions 
are affribufed fo fhe frilepfon channel and removed from fhe lefralepfon channel. The dilepfon channels 
are nol explicilly required to be orlhogonal wifh fhe fefralepfon channel, buf fhe overlap is found to be 
negligibly small in simulafed samples and non-exislenl in dafa. 

Background evenfs confaining well-identified prompf leptons are modelled by simulation. The norm- 
alisalions used for fhe backgrounds in Ihis calegory are laken from dafa conlrol regions if fhe resulting 
normalisalion uncerfainly is lower fhan fhaf from fhe Iheorelical prediction. The yields in fhe dafa conlrol 
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regions are extrapolated to the signal regions using the simulation. Background sources involving one or 
more incorrectly identified lepton, e.g. instrumental backgrounds, are modelled using data events from 
control regions, except in the opposite-sign dilepton channel, where this background is very small. 

The following sections describe additional selection requirements and the background evaluation in each 
of the four channels. 


5.1 Opposite-sign dilepton channel 

In the opposite-sign dilepton channel, events are required to have at least three jets, one or two of which 
are f^-tagged. Two orthogonal selections are defined fo separafe ttW and ttZ final sfafes. The firsl (2f’-noZ) 
selecfs differenf-llavour lepfon (DF) evenfs wifh fhe scalar sum of fhe px of lepfons and Jefs, Hj, above 
130 GeV, and same-flavour lepfon (SF) evenfs fhat are nof compatible wifh Z boson or low-mass res¬ 
onance producfion, by requiring \ma - mz\ > 10 GeV, ma > 15 GeV and > 40 GeV. The 2GZ 
selection contains SF events within the mass window Inia - mz\ < 10 GeV. In both selections an addi¬ 
tional requirement on the average distance between two jets, calculated using all possible jet pairs in the 
event, > 0.75, is applied to remove the low-dijet-mass region where the Alpgen-i-Pythia simulation 
does not provide a good description of the Z boson background [67] . 

For 27’-noZ events, the ttV signal contribution originates mainly from the tt dilepton final sfafe accom¬ 
panied by a hadronic W/Z boson decay and from fhe tt single-lepfon final sfafe wifh a lepfonic W boson 
decay. For 27’-Z evenfs, fhe confribufion of tiW producfion is negligible while fhe ttZ confribufion comes 
from fhe fully hadronic ti final sfafe wifh a lepfonic Z boson decay. 

Affer evenf selecfion fhe dominanf backgrounds are tt and Z producfion in l£-noZ and 2C-Z respecf- 
ively, and fhe exfracfion of fhe signal relies on discriminating if from fhese backgrounds, based on well- 
modelled evenf kinemafics. To improve fhe modelling of fhe tt background, fhe simulafed tt evenfs are 
reweighfed fo accounf for fhe observed differences in fhe fop quark px and fhe tt sysfem px befween 
dafa and Powheg-i-Pythia simulafion in measuremenfs of differential cross secfions at ^[s -1 TeV [68]. 
To improve Z background modelling, fhe simulafion is reweighfed fo accounf for fhe difference in fhe 
Z Px specfrum befween dafa and simulafion [67], and fhe ZQQ (Q - b,c) componenf of fhe Z back¬ 
ground is adjusfed fo mafch dafa in a ZQQ-dominafed confrol region wifh af leasf one f7-fagged jef. Small 
background confribufions arise from single-fop-quark Wt channel producfion, diboson (WW, VTZ, ZZ) 
processes, fhe associafed producfion of a Higgs boson and a tt pair, fhe associafed producfion of a WW 
and a tt pair, and fhe associafed producfion of a single fop quark and a Z boson. All of fhese backgrounds 
are defermined from simulafion. 

In fhe 2^-noZ region, W boson, tt (wifh a single lepfon in fhe final sfafe) and t- and ^-channel single-fop- 
quark producfion processes can safisfy fhe selecfion requiremenfs due fo fake lepfons. These backgrounds 
are a small fraclion of fhe fofal esfimafed background, and fheir yields are esfimafed using simulafion and 
cross-checked wifh a data-driven technique based on the selection of a same-sign lepton pair. 

Events are categorised according to the number of jets and the number of f^-tagged jets. In the 27’-noZ 
selection, events with one or two f^-tagged jets are separated into three exclusive regions according to the 
jet multiplicity, with three (27’-noZ-3j), four (27’-noZ-4j), and five or more (27’-noZ-5j) jefs. In fhe 2£-Z 
selecfion, events with exactly two f^-tagged jets are separated into three regions according to the same 
scheme: 27’-Z-3j, 27’-Z-4j and 27’-Z-5j. 


A neural network (NN) discriminant built using the NeuroBayes [69] package is used to separate the 
combined ttW and ttZ signal from the background in the signal-rich regions 2^-noZ-4j, 2^-noZ-5j and 
2AZ-5j. The other regions considered in the opposite-sign dilepton channel have lower sensitivity and 
are used as control regions; event counting is used in the 2f’-noZ-3j region, while the scalar sum of the jet 
transverse momenta is used as a discriminant in the 2^-Z-3j and 2AZ-4j regions. The inclusion of 
these highly populated control regions, enriched in tt or Z backgrounds, in the fit used to extract the ttV 
signals, strongly constrains the normalisation uncertainties of these backgrounds. This in turn improves 
the background predictions in the signal-rich regions. The signal and control regions are summarised in 
Table 2. 


Region 

Targeting 

Sample fraction [%] 

2AnoZ-4j 

ttW and ttZ 

0.68 

2AnoZ-5j 


1.2 

2AZ-5j 

ttZ 

3.3 

27’-noZ-3j 

tt 

92 

2AZ-3j 

z 

70 

2AZ-4j 


66 


Table 2; Signal and control regions of the opposite-sign dilepton channel, together with the processes targeted and 
the expected fraction of the sample represented by the targeted process. 

The set of variables used as input to the NN discriminant is chosen separately for each signal region, 
based on the ranking procedure implemented in the NeuroBayes package which takes into account the 
statistical separation power of the variables and the correlations between them. All variables used for 
the NN training are required to show good agreement between data and background expectation in the 
control regions. Seven variables are selected in each signal region. The list of selected variables and their 
ranking is shown in Table 3. 

Figure 2 illustrates the discrimination between the ttV signal and background provided by the NN dis¬ 
criminants. Since in the 2f’-noZ region the contributions from both ttW and ttZ production are comparable 
in size and have similar kinematics, they result in a similar NN discriminant shape and are thus fitted to¬ 
gether. In the 2£-X region, the tfW contribution is negligible, and thus the NN discriminant shape is driven 
by the ttZ signal. 

The expected sample compositions in each of the three signal and three control regions are summarised 
in Table 4 along with the number of events observed in data. The distributions of discriminants in the 
control regions are shown in Figure 3. The data and simulation agree within the expected uncertainties. 


5.2 Same-sign dilepton channel 

The same-sign dilepton channel targets the ttW process. Events are required to have E™®® > 40 GeV, 
//t > 240 GeV and to contain at least two f7-tagged jets. The same-sign dilepton channel is divided into 
three orthogonal regions based on the flavour combination of fhe lepton pair: 2e-SS, efi-SS, and 2/r-SS. In 
fhe 2e-SS region, an addifional requiremenf on fhe dilepfon mass removing evenfs wifh 75 GeV < nige < 
105 GeV is imposed fo reduce fhe contamination by Z —> ee events where the charge of one electron is 
misidentified. A similar requirement is not imposed on the efi-SS or 2fi-SS regions, since the probability 
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NN rank 


Variable 

Definition 

27’-noZ-4j 

2GnoZ-5j 

2GZ-5j 


Invariant mass of the two highest pj untagged 
jets in events with exactly two fj-tags, or of the 
two highest pj untagged jets, excluding the jet 
with the second highest b-tag weight, in events 
with exactly one b-tag 

1st 

7th 


Centralityjgj 

Sum of pj divided by sum of E for all jets 

2nd 

1st 

6th 

Hi 

2nd Fox-Wolfram moment [70] 

3rd 

2nd 

- 

min AR 

jj 

Invariant mass of the combination of the two 
jets with the smallest A/? 

4th 

6th 

- 

max m™" ^ 

Larger of the invariant masses of the two 
(lepton, (j-tagged jet) pairs, which are built 
based on the minimum A7?(f, b) for each lepton 

5th 

5th 



Jet3 

rj 

Third-jet px 

6th 

- 

- 

Jet4 

rj 

Fourth-jet px 

- 

3rd 

- 

A/?"e 

Average A/? for all jet pairs 

7th 

- 

- 

\mjj-mv\<30 

jet 

Number of jet pairs with mass within a 30 GeV 
window around 85 GeV 

- 

4th 

2nd 

■''^40 

Number of jets with px > 40 GeV 

- 

- 

1st 

max pT 
^bb 

Invariant mass of the combination of two b- 
tagged jets with the largest vector sum px 

- 

- 

3rd 

^tit2 

A/? between the two leptons 

- 

- 

4th 

b] 

Invariant mass of the combination of the two 
jets with the largest vector sum px; one jet must 
be (j-tagged 



5th 


2nd Fox-Wolfram moment built from only jets 

- 

- 

7th 


Table 3: Definitions and rankings of the variables considered in each of the regions where a NN is used in the 
opposite-sign dilepton channel. 
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Figure 2: The NN output distributions for the three signal regions in the opposite-sign dilepton channel, before the 
fit to data. The distributions are shown in the (a) 2f’-noZ-4j, (b) 2^-noZ-5j and (c) 2f’-Z-5j regions. The orange 
dashed lines show the tIV signal normalised to the background yield. “Rare SM” comprises the diboson, single-top, 
tZ, WtZ, ttH processes and the fake lepton background. The hatched area corresponds to the total uncertainty on 
the predicted yields. The “Data/SM” plots show the ratio of the data to the total Standard Model expectation. 



Figure 3: Control region distributions in the opposite-sign dilepton channel, before the ht to data. The distributions 
are shown in the (a) 2^-noZ-3j, (b) 2f’-Z-3j and (c) 2^-Z-4j regions. For the 2^-noZ-3j region only the event count is 
used in the fit. The hatched area corresponds to the total uncertainty on the predicted yields. “Rare SM” comprises 
the diboson, single-top, fZ, WtZ, ttH processes and the fake lepton background. The “Data/SM” plots show the 
ratio of the data to the total expected Standard Model expectation. 
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Region 

t + X 

Bosons 

Fake leptons 
charge mislD 

Total expected 
background 

ttW 

ttZ 

Data 

2AnoZ-3j* 

20800 + 2600 

600 + 200 

160 + 80 

21600 + 2700 

42.0 + 2.8 

23.2+1.5 

22585 

2f’-noZ-4j 

8200 + 1400 

240 + 90 

80 + 40 

8600 +1400 

36.6+1.8 

22.4+1.1 

8909 

2f’-noZ-5j 

3700 + 850 

100 + 40 

47 + 23 

3810 + 870 

24.9 + 2.2 

22.4 + 2.0 

3901 

2AZ-3j* 

800 + 140 

1960 + 880 

4.1 + 2.1 

2760 + 890 

1.24 + 0.13 

3.71 + 0.38 

2806 

2AZ-4j* 

330 + 70 

740 + 390 

2.2+ 1.1 

1100 + 400 

1.31+0.11 

7.21 + 0.58 

1031 

2CZ-5j 

170 + 40 

340 + 200 

1.4+ 0.7 

510 + 210 

0.89 + 0.07 

17.7+1.4 

471 

2e-SS 

0.66 + 0.13 

0.17 + 0.10 

8.9+ 2.4 

9.8+ 2.6 

2.97 + 0.30 

0.93 + 0.23 

16 

gyU-SS 

1.9 + 0.35 

0.39 + 0.28 

14.1+4.5 

16.4 + 5.1 

8.67 + 0.76 

2.16 + 0.51 

34 

2//-SS 

0.94 + 0.17 

0.25 + 0.14 

0.93 + 0.55 

2.12 + 0.86 

4.79 + 0.40 

1.12 + 0.27 

13 

3AZ-0b3j* 

1.11 + 0.32 

67+16 

15.2 + 6.0 

83+15 

0.05 + 0.03 

1.86 + 0.47 

86 

3AZ-lb4j 

1.58 + 0.42 

3.8+ 1.3 

2.4+ 1.1 

7.8+1.6 

0.14 + 0.05 

7.1 + 1.6 

8 

3AZ-2b3j 

1.29 + 0.34 

0.68 + 0.33 

0.19 + 0.13 

2.16 + 0.42 

0.21+0.07 

2.76 + 0.69 

3 

3AZ-2b4j 

1.00 + 0.29 

0.48 + 0.24 

0.42 + 0.37 

1.93 + 0.49 

0.14 + 0.07 

6.6+ 1.6 

11 

3f’-noZ-2b 

1.06 + 0.25 

0.27 + 0.17 

1.31 + 0.90 

2.7+ 0.9 

3.7+ 0.9 

1.23 + 0.32 

6 

4CDF-0b 

0.06 + 0.01 

0.11 + 0.04 

0.03 + 0.17 

0.21+0.22 

- 

0.28 + 0.01 

2 

4CDF-lb 

0.22 + 0.03 

0.05 + 0.03 

0.13 + 0.22 

0.39 + 0.27 

- 

1.05 + 0.03 

1 

4CDF-2b 

0.11 + 0.02 

<0.01 

0.11 + 0.19 

0.22 + 0.21 

- 

0.64 + 0.02 

1 

4CZZ* 

0.01 + 0.00 

134.2+1.2 

0.27 + 0.18 

134.5+1.3 

- 

0.07 + 0.01 

158 

4ASF-lb 

0.16 + 0.02 

0.29 + 0.06 

0.14 + 0.19 

0.61+0.27 

- 

0.91 + 0.02 

2 

4ASF-2b 

0.08 + 0.01 

0.09 + 0.03 

0.04 + 0.18 

0.21+0.23 

- 

0.64 + 0.02 

1 


Table 4: Expected event yields for signal and backgrounds, and the observed data in all signal and control regions 
(marked with an asterisk) used in the fit to extract the ttW and ttZ cross sections. The quoted uncertainties on 
expected event yields represent systematic uncertainties including MC statistical uncertainties. The tt, single-top, 
tZ, WtZ, tlH and tirf processes are denoted t + X. The Z, WW, fkZ, ZZ, ttWW and processes are denoted 

‘Bosons’. 


for the muon charge to be misidentified is found to be negligible, and Z-i-jets is not a dominant background 
in the efi region. 

Signal events from the ttW process are produced when the associated W boson decays leptonically and 
the tt system decays in the f’-i-jets channel. 

A smaller contribution from ttZ comes from a leptonic decay of the Z boson where one lepton is not 
reconstructed, together with a leptonic decay of one of the two W bosons coming from the top quark 
decays. 

The main backgrounds vary depending on the lepton flavour: events containing a lepton with mis¬ 
identified charge are dominant in the 2e-SS region and prevalent in the e/r-SS region, whereas events 
with a fake lepton contribute significantly in all regions, but are dominant in the 2jU-SS region. Back¬ 
grounds from the production of prompt leptons with correctly identified charge come primarily from WZ 
production, but these are small compared to the instrumental backgrounds. 

Processes featuring an opposite-sign lepton pair, like tt and Z boson production, can enter this channel 
through the misidentification of the electron charge. Charge misidentification rates, parameterised in 
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Pt and \r]\ of the electrons, are measured in a control region containing events with two electrons with 
75 < niee < 105 GeV, which is divided into same-sign and opposite-sign subregions with non-Z ^ ee 
backgrounds subtracted. A likelihood function is constructed relating the number of observed events in 
the two subregions with the probability for an electron falling in a given (pT^ I 7 I) bin to be reconstructed 
with the wrong charge, and maximised to obtain the charge misidentification rates. 

A template is then constructed using opposite-sign data with event selection identical to that used in the 
signal region except for the requirement on the charge of the leptons. A weight given by 

e\ + S2- 2SiE2 

'" = 1—i—I-^^ 

1 - (£1 -I- £2 - 2£i£2) 

is applied to each event and used to construct the template, where the charge misidentification rates for 
the two leptons are £ 12 - These are set to zero in the case of muons. 

To estimate the background from fake leptons in the same-sign dilepton channel a set of scale factors are 
measured. The scale factors are defined as f = Nj/Nj^, fhe ratio of fhe number of observed fighf lepfons, 
i.e. lepfons safisfying all selecfion criferia, fo fhe number of loose lepfons. Loose lepfons differ from fighf 
lepfons in fhaf fhey are required fo fail isolafion requiremenfs; loose muons additionally have relaxed 
selection criferia, requiring |zol < 2 mm wifh no requiremenf on do. The scale facfors are measured in a 
confrol region (2^-SS-CR) confaining fwo same-sign lepfons (vetoing evenfs wifh a fhird lepfon, as is done 
in fhe signal region), af leasf one Z?-fagged jef, and Hj < 240 GeV. The missing fransverse momenfum 
disfribufions in fhis confrol region are shown in Figure 4. A femplafe for fake lepfon backgrounds is 
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Figure 4: Distributions of F™'*® for events in a same-sign dilepton control region with Hj < 240 GeV and at least 
one f7-tagged jet for the different lepton flavour combinations, (a) ee, (b) ep and (c) fufu. “Rare SM” contains small 
background contributions mainly consisting of the WW and WZ processes. The predictions are shown before fitting 
to data in the control region. The instrumental backgrounds, including fake leptons and leptons with misidentified 
charge are predicted using data-driven methods. The hatched area corresponds to the total uncertainty on the 
predicted yields. The last bin in each histogram includes the overflow. 

constructed in the control region from the loose lepton sample, using the expression 
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where A^tl is the number of events in which the first lepton is a tight lepton and the second is loose, and 
^LT, ^LL are defined in a similar fashion. The fi^ are the scale factors for the first and second leptons. 
A subtracted term, shown in the brackets labelled “MC, prompt” is included to remove contamination 
from prompt lepton production in the loose lepton sample. This subtraction accounts for about 2-3% of 
the total estimate. The background template for fake leptons is fitted to data in the control region in bins 
of lepton pt to obtain the scale factors ft, which are measured separately for electrons and muons and 
also binned in lepton pj. With the measured scale factors, the background templates are produced in the 
signal region according to Eq. (2). 

Potential overlap between the estimates of charge misidentification and fake leptons is taken into account 
with an additional subtraction step. A background template for fake leptons is produced for opposite-sign 
events (using a selection and binning otherwise the same as the signal region), and the charge misidentific¬ 
ation rates are applied to this template to obtain a representation of the overlap of these two backgrounds. 
This new template is then subtracted from the total estimate of the fake lepton background. This subtrac¬ 
tion represents about 2%-5% of the total yield. 

To improve the separation of the ttW signal from backgrounds, events in the ep-SS and 2p-SS regions 
are further divided into four bins based on jet multiplicity and missing transverse momentum. Events are 
classified as low-Ajets (Ajets - 2 or 3) or high-Ajets (Ajets > 4) and low-E™'*® (40 < < 80 GeV) or 

high-E™*® (E™®* > 80 GeV). No further event classification is used in the 2c-SS region. The expected 
and observed contributions in each of the three dilepton flavour regions are summarised in Table 4 and 
plotted in Eigure 5. 



Figure 5; Event yields in the same-sign dilepton signal regions according to the binning used in the final likelihood 
fit in the (a) 2e-SS (b) ep-SS and (c) 2p-SS regions. The disttibutions are shown before the fit. The bins labelled 
“Low-E™'**” correspond to 6 (40,80) GeV, and those labelled “High-E™"’*” correspond to E™*** > 80 GeV. 
“Rare SM” contains small background contributions mainly consisting of the ttH and ITZ processes. Instrumental 
backgrounds, including fake leptons and leptons with misidentified charge are predicted using data-driven methods. 
The hatched area corresponds to the total uncertainty on the predicted yields. 


5.3 Trilepton channel 

In the trilepton channel two preselections are considered, referred to as and 3^-noZ. The 3^-Z region 
targets the ttZ process, while 3f’-noZ aims at measuring the ttW process. In the region 3GZ, at least one 
pair of leptons is required to have the opposite sign and same flavour (OSSE) and have an invariant 
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mass within lOGeV of the Z boson mass. Region 3^-noZ contains the remaining trilepton events with a 
requirement that the leptons must not all have the same sign. 

The trilepton channel signal regions are determined as follows. First, the preselected samples are split into 
categories according to the jet multiplicity and the number of Zj-tagged jets. The categories with similar 
predicted signal-to-background ratio (S/B) and systematic uncertainties are grouped together. The final 
selection in each group is optimised for maximal expected significance, including bofh fhe sfafisfical and 
sysfemafic uncerfainfies, using requiremenfs on and lepfon pj. If is found fhaf opfimal significance 
is obfained wifhouf a requiremenf on 

Four signal regions are defined as a resulf of fhe grouping and opfimisafion: 3GZ-lb4j, 3GZ-2b3j, 
3f’-Z-2b4j and 3GnoZ-2b. In fhe 3f’-Z-lb4j region, af leasf four jefs are required, exacfly one of which is 
b-tagged. In fhe 3GZ-2b3j region, exacfly fhree jefs wifh af leasf fwo f^-tagged jefs are required. In fhe 
3GZ-2b4j region af leasf four jefs are required, of which af leasf fwo jefs are b-tagged. In fhe 3GnoZ-2b 
region af leasf fwo and af mosf fhree jefs are required, of which af leasf fwo jefs are f^-fagged. For evenfs 
in which fhe fhird leading lepfon is an elecfron, fhe minimum pj requiremenf on fhe fhird lepfon is raised 
fo 20GeV in fhe 3GZ-lb4j, 3GZ-2b3j and 3GZ-2b4j regions, and fo 25 GeV in fhe 3f’-noZ-2b region. 

The 3{-Z preselecfion is dominafed by fTZ evenfs, wifh a significanf confribufion from evenfs wifh fake 
leptons. To consfrain fhe fTZ background, a confrol region called 3^-Z-0b3j is defined and included in fhe 
fif. In fhis region, fhe presence of exacfly fhree jefs, wifh exacfly zero b-tags, is required in addition to fhe 
requiremenfs of fhe 3£-Z preselecfion. The normalisafion correction for fhe fTZ background wifh respecf 
fo fhe Sfandard Model expecfafions is obfained from fhe fif and found fo be 0.98 ± 0.20. The quofed 
uncerfainfy includes bofh fhe sfafisfical and sysfemafic componenfs. The modelling of WZ production 
in association wifh heavy-flavour jefs is furlher validated in a confrol region 3GZ-lb-CR, defined by 
requiring fhe presence of one fo fhree jefs, exacfly one of which is b-tagged. 

The fake lepfon background is esfimafed by using fhe so-called mafrix mefhod [71], which makes use of 
an orthogonal confrol region in which lepfon isolation and elecfron idenfificafion criferia are relaxed. The 
efficiencies for real and fake lepfons used in fhe mafrix mefhod are measured in evenfs confaining fwo 
leptons and one f^-fagged jef. To validafe fhe esfimafe of fhe background confaining fake lepfons, a confrol 
region 3f’-noZ-lb-CR is defined by requiring exacfly one jef fo be fj-fagged in addifion fo fhe requiremenfs 
of fhe 3f’-noZ region. Figure 6 shows disfribufions of and fhird-lepfon pj in fhe 3^-noZ-lb-CR and 
3f’-Z-lb-CR regions, respecfively. The level of agreemenf befween dafa and expecfafion is good. 

The signal and confrol regions of fhe frilepfon channel used in fhe fif are summarised in Table 5. The 
expecfed and observed yields in fhe signal and confrol regions are shown in Table 4. Even! yields sum¬ 
marising fhe signal regions wifh differenl lepton flavour combinations (3f’-Z-SR) and fhe disfribufion of 
fhe minimum invarianf mass of jef friplefs (minimum injjj) for evenfs in fhe 3(’-Z-2b4j region are shown in 
Figure 7. Considering fhe four leading jefs, fhe momenfum vecfor sum of fhe minimum invarianf jef friplef 
mass is found fo give a powerful esfimafe of fhe hadronically decaying fop direcfion. Good agreemenf 
befween dafa and expecfafion is observed. 


5.4 Tetralepton channel 

The fefralepfon channel fargefs fhe ttZ process for fhe case where bofh W bosons resulfing from top quark 
decays and fhe Z boson decay lepfonically, and uses an even! counfing approach in five signal regions. 
Evenfs wifh fwo pairs of opposife-sign lepfons are selecfed, among which af leasf one pair is same flavour. 
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Figure 6; Distributions of (a) in the Sf-noZ-lb region and (b) third-lepton pj in the 3^-Z-lb region. “Rare 
SM” contains small background contributions consisting of the WWW, WWZ, H —> ZZ and ttWW processes. The 
hatched area corresponds to the total uncertainty on the predicted yields. The distributions are shown before the fit. 
The last bin in each histogram includes the overflow. 


Region 

Targeting 

Sample fraction [%] 

3GZ-lb4j 


M 

3^-Z-2b3j 

ttZ 

54 

3^-Z-2b4j 


76 

3GnoZ-2b 

ttW 

48 

3^-Z-0b3j 

WZ 

68 


Table 5; Signal and control regions of the trilepton channel used in the fit, together with the processes targeted and 
the expected fraction of the sample represented by the targeted process. 


The OSSF lepton pair with reconstructed invariant mass closest to mz is attributed to the Z boson decay 
and denoted in the following as Z\. The two remaining leptons are used to define Z 2 . The signal regions 
are defined according fo fhe relative flavour of fhe fwo remaining lepfons, differenf or same flavour, and 
fhe number of f^-tagged jefs: zero, one, or af leasf fwo {Ob, lb, 2b). The signal regions are fhus 4^-DF-Ob, 
4f’-DF-lb, 4f’-DF-2b, 4^-SF-lb and 4^-SF-2b. The ZZ background mosfly afFecls fhe SF regions and 
Iherefore evenfs wifh a Z 2 SF lepton pair and no f^-lagged jefs are discarded. 

Furfher requiremenfs are applied in each signal region such fhaf fhe expecfed sfafisfical uncerfainfy on 
fhe measured ttZ signal cross secfion is minimised. Evenfs in fhe 4^-SF-lb region are rejecfed if fhey are 
compatible wifh a ZZ even!, i.e. if E™®® < 80(40) GeV for mz^ inside (oufside) a 10 GeV region cenfred 
af fhe Z boson mass. This requirement on is relaxed by 40 GeV for the 4GSF-2b region. The 
impact of events with fake leptons decreases with the number of reconstructed (^-tagged jets. To suppress 
these backgrounds, additional requirements on the scalar sum of the transverse momenta of the third and 
fourth leptons (pt 34 ) are imposed in the lower fj-tag multiplicity regions. In the TGSF-lb, 4GDF-lb and 
TGDF-Ob regions events are required to satisfy pt 34 > 25 GeV, pt 34 > 35 GeV and pt 34 > 45 GeV, 
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Figure 7: (a) Event yields in the trilepton channel summarising the signal regions with different lepton flavour 
combinations and (b) the minimum three-jet invariant mass for events in the 37’-Z-2b4j signal region. “Rare SM” 
contains small background contributions consisting of the WWW, WWZ, H ZZ and ttWW processes. The 
distributions are shown before the fit. The hatched area corresponds to the total uncertainty on the predicted yields. 
The last bin in (b) includes the overflow. 


respectively. In the 4^-DF-Ob region the requirement on the fourth lepton is raised to pj > 10 GeV 
and at least two jets must be reconstructed in the event. In all regions, the invariant mass of any two 
reconstructed OS leptons is required to be larger than 10 GeV. The definitions of the signal regions are 
summarised in Table 6. 


Region 

Z 2 leptons 

Pta 

PT 34 

|m« -mzjl 

^miss 

^Yiets 

^&-jets 

4f’-DF-0b 


> 10 GeV 

> 45 GeV 

- 

- 

> 2 

0 

4GDF-lb 


> 7 GeV 

> 35 GeV 

- 

- 

- 

1 

4f’-DF-2b 


> 7 GeV 

- 

- 

- 

- 

> 2 

4GSF-lb 


> 7 GeV 

> 25 GeV ■ 

f > 10 GeV 
< 10 GeV 

> 40 GeV 'I 

> 80 GeV , 

- 

1 

4GSF-2b 


> 7 GeV 

- i 

f > 10 GeV 
< 10 GeV 

> 40 GeV , 

- 

> 2 


Table 6: Definitions of the five signal regions in the tetralepton channel. 


The ZZ background is large in the tetralepton channel, and therefore a control region, A£-7X, is defined fo 
consfrain fhe ZZ normalisafion in fhe SF region, and is included in fhe fil. Bofh lepfon pairs are required 
fo salisfy |niz, 2 ~ ^z\ < 10 GeV, and evenfs are refained if F"™*® < 50 GeV. The filled normalisafion 
correction wilh respecf fo fhe Sfandard Model expeclalion is 1.16 + 0.12. The quoted uncerlainlies include 
bofh fhe slalisfical and syslemalic componenls. The number of jels and fi-lagged jels in fhe df-ZZ region 
are shown in Figure 8. Dala disfribulions agree wifh expecfalions from simulalion. 

The telralepton channel backgrounds wilh al leasl one fake lepfon are esfimaled using simulalion, where 
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Figure 8: Distributions of (a) the number of jets and (b) the number of h-tagged jets in the ZZ control region of 
the tetralepton channel 4E-ZZ. The hatched area corresponds to the total uncertainty on the predicted yields. The 
distributions are shown before the fit. “Rare SM” contains small background contributions mainly consisting of 
the WtZ and tJH processes. The “Data/SM” plots show the ratio of the data events to the total Standard Model 
expectation. 


the prediction is corrected with a constant factor to improve agreement in control regions with enhanced 
fraction of single fake leptons. By probing the fake muon and electron background in a Z+fake-lepton- 
candidate control region and in a tt+fake-lepton-candidate region, two equations per lepton flavour / can 
be constructed and the correction factors for these two processes can be determined: c'' = 1.23 ±0.13, 
c” = 1.25 ± 0.09, cf = 1.35 ± 0.05, and = 1.61 ± 0.05. The quoted uncertainties include systematic 
effects. The control regions are required to contain three leptons and are either Z- or ri-like. In the first 
case an OSSF lepton pair is required together with < 30 GeV. In addition, the transverse mass of 
the non-Zi lepton i is required to satisfy mj < 35 GeV, where niT is defined as 

mj - (3) 

The non-Zi lepton is then used as the fake lepton candidate. In the second case an OS lepton pair is 
required together with at least one jet with pj > 30 GeV, and events with an OSSF lepton pair are 
rejected. The lowest-px same-sign lepton is then used as the candidate. The background from events 
with two fake leptons is evaluated from simulation with relaxed requirements and extrapolated in several 
steps into the signal region. The total background yield and its uncertainty are dominated by the estimate 
extracted from simulation of trilepton events with only one additional fake lepton. 

The expected sample composition of the six tetralepton regions is summarised in Table 4 along with the 
number of events observed in data. Seven events are observed in the five signal regions (4L-SR). Figure 9 
shows good agreement between data and expectation for the distributions of the number of jets, number 
of Z?-tagged jets, as well as the invariant masses of the two pairs of leptons. 
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Figure 9; Distributions of (a) number of jets, (b) number of b-tagged jets, invariant mass of the (c) Zi and (d) Z 2 
dilepton pair for the tetralepton signal region selection. The distributions are shown before the fit. The distributions 
of the seven observed events are compared to expectation. The hatched area corresponds to the total uncertainty 
on the predicted yields. “Rare SM” contains small background contributions mainly consisting of the WtZ and tIH 
processes. 
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6 Systematic uncertainties 


Several sources of systematic uncertainty are considered that can affect the normalisation of signal and 
background in each channel and/or the shape of the discriminant distributions in the opposite-sign dilepton 
channel. 

The luminosity estimate has an uncertainty of 2.8%, determined using beam-separation scans [72]. This 
systematic uncertainty is assigned to all background contributions obtained from MC simulation. 

6.1 Uncertainties on reconstructed objects 

Uncertainties associated with the lepton selection arise from the imperfect knowledge of the trigger, 
reconstruction, identification and isolation efficiencies, and lepton momentum scale and resolution. The 
uncertainty on the electron identification efficiency is the largest systematic uncertainty in the trilepton 
channel and among the most important ones in the tetralepton channel. 

Uncertainties associated with the jet selection arise from the jet energy scale (JES), JVF requirement, jet 
energy resolution (JER) and jet reconstruction efficiency. The JES and its uncertainty are derived com¬ 
bining information from test-beam data, collision data and simulation [63]. JES uncertainty components 
arising from the in-situ calibration and the jet flavour composition are among the dominant uncertainties 
in the opposite-sign dilepton, same-sign dilepton and trilepton channels. The uncertainties in the JER and 
JVE have a significant effect at low jet pT- The JER uncertainty is the second largest uncertainty in the 
trilepton channel while the JVE uncertainty is not negligible in the opposite-sign dilepton, trilepton and 
tetralepton channels. 

The efficiency of the flavour tagging algorithm is measured for each jet flavour using control samples in 
data and in simulation. From these measurements, correction factors are defined to correct the tagging 
rates in the simulation. In the case of fj-jets, correction factors and their uncertainties are estimated based 
on observed and simulated fj-tagging rates in tt dilepton events [65]. In the case of c-jets, they are derived 
based on jets with identified D* mesons [73]. In both cases the correction factors are parameterised as a 
function of jet pj- In the case of light-flavour jets, correction factors are derived using dijet events, and 
are parameterised as a function of jet pT and rj [73]. Sources of uncertainty affecting the b- and c-tagging 
efficiencies are considered as a function of jet px, including bin-to-bin correlations [65]. An additional 
uncertainty is assigned to account for the extrapolation of the fj-tagging efficiency measurement from the 
Px region used to determine the scale factors to regions with higher px. For the light-jet tagging efficiency 
the dependence of the uncertainty on the jet px and rj is considered. These systematic uncertainties are 
taken as uncorrelated between fj-jets, c-jets, and light-flavour jets. 

The treatment of the uncertainties on reconstructed objects is common to all four channels, and thus these 
are considered as correlated among different regions. 


6.2 Uncertainties on signal modelling 

To assess the factorisation and renormalisation scale uncertainties on ttV modelling, the scales are varied 
up and down by a factor of two in a correlated manner. Radiation uncertainties are assessed by simul¬ 
taneously varying the scale of the momentum transfer Q in the running strong coupling as(Q^) in the 
matrix-element calculation and in the Pythia parton shower, up or down by a factor of two [74]. 
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In addition, the jet pj matching threshold and the amount of radiation in the parton shower are inde¬ 
pendently varied up and down by a factor of two. The dominant systematic uncertainty comes from the 
variation of Q in as(Q^) in the matrix element calculation and in the Pythia parton shower. This variation 
has a significant effect on the distribution of the number of jets in tiV events. 

Systematic uncertainties due to the choice of PDF are evaluated using the uncertainty sets of the CTIO 
NLO, MSTW2008 68% confidence level (CL) NLO and NNPDF 2.3 NLO [75] PDFs following fhe 
PDF4LHC recommendafions [76]. 

The uncerfainfies on fhe ttV modelling are among fhe dominanf ones in fhe fefralepfon channel buf fhey 
have a negligible impacf in all ofher channels. Signal modelling and PDF uncerfainfies are freafed as 
correlafed among channels. 


6.3 Uncertainties on background modelling 

Uncerfainfies on fhe background modelling differ significanfly among fhe channels due fo large differ¬ 
ences in fhe background composifion. 


Z boson background: This dominafes in fhe 2€-Z regions of fhe opposife-sign dilepfon channel. Four 
sources of uncerfainfy are considered: fhose associafed wifh fhe cross section, fhe Z boson pj correcfion, 
fhe scale choice for parfon emission, and fhe choice of generafor, evaluafed by comparing fhe nominal 
Alpgen sample fo a Sherpa sample generafed using Sherpa 1.4.1 wifh up fo fhree addifional parfons in fhe 
LO mafrix elemenf and fhe CTIO PDF sef. 


tt background: This dominafes in fhe 2^-noZ regions of fhe opposife-sign dilepfon channel. A number 
of sysfemafic uncerfainfies affecting fhe modelling of fhe ti process are considered in Ibis channel: fhose 
due fo fhe uncerfainfy on fhe cross section which amounf fo +5%l-6%, due fo fhe choice of parfon shower 
and hadronisafion model (evaluafed by comparing evenfs produced by Powheg inferfaced wifh Pythia or 
Herwig [77]), due fo fhe choice of generafor (evaluafed by comparing a sample generafed using Madgraph 
inferfaced wifh Pythia fo fhe defaulf tt sample), and due fo fhe reweighfing procedure applied fo correcf 
fhe ttMC modelling. An addifional 50% normalisafion uncerfainly is assigned fo tt -i-heavy-flavour (HF) 
jefs production fo accounf for limifed knowledge of Ibis process. 


Single-top background: This is small and affecls only fhe opposife-sign dilepfon channel. An un¬ 
certainly of 6.8% is assigned fo fhe cross secfion for single-fop producfion [42], corresponding fo fhe 
Iheorelical uncerfainly on Wt producfion, fhe only process conlribufing fo fhis final slafe. An addifional 
conlribufion arises from fhe comparison of predictions using differenl schemes fo accounf for inferference 
befween Wt and tt. 
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Diboson background: In the trilepton and same-sign dilepton channels the diboson background is 
dominated by IVZ production, while in the tetralepton channel ZZ production is dominant. In the 
opposite-sign dilepton channel the diboson background includes WW, 'WZ and ZZ production and the 
uncertainties are assigned to the sum of these processes. 

In the trilepton and same-sign dilepton channels, the normalisation of the 'WZ background is treated as a 
floating parameter in the fit used to extract the ttV signal. The uncertainty on the extrapolation of the VTZ 
background estimate from the control region to signal regions with specific jet and Z^-tag multiplicities 
is evaluated by comparing the nominal Sherpa sample to the prediction of Powheg, as well as by using 
variations of the simulation parameters. The uncertainty amounts to 20%-35%. 

In the tetralepton, trilepton and same-sign dilepton channels the normalisation of the ZZ background is 
treated as a floating parameter in the fit used to extract the ttV signal. In the tetralepton channel, several 
uncertainties on the ZZ background estimate are considered. They arise from the extrapolation from the 
A{-ZZ control region (corresponding to on-shell ZZ production) to the signal region (with off-shell ZZ 
background) and from the extrapolation from the control region without jets to the signal region with at 
least one jet. Using data-driven techniques, these uncertainties are found to be 30% and 20%, respectively. 
An additional uncertainty of 10%-30% is assigned to the normalisation of the heavy-flavour content of 
the ZZ background based on a data-to-simulation comparison of events with one Z boson and additional 
jets, and cross-checked with a comparison between different ZZ simulations. 

In the opposite-sign dilepton channel, in which the diboson background is small, a 20% uncertainty is 
assigned to the VTZ and ZZ background normalisation. This is estimated from the level of agreement 
between data and prediction in the 3^-Z-0b3j control region. 


ttH background: An uncertainty of 12% is assigned to the ttH production cross section [47] in all 
channels. Additional uncertainties that affect ttH kinematics are assigned in the opposite-sign dilepton 
channel and are negligible for the other channels. These uncertainties come from the choice of factorisa¬ 
tion and renormalisation scales, and the functional form of the scale in ttH samples. 


tZ and WtZ background: In the opposite-sign dilepton and trilepton channels, tZ and WtZ back¬ 
grounds are summed and an uncertainty of 20% is assigned to their cross section. An additional uncer¬ 
tainty on the shape of this background is considered in the trilepton channel for which this background 
is important. The shape uncertainty is evaluated by varying the factorisation and renormalisation scales 
and as in simulation up and down by a factor of two with respect to the nominal value, in a correlated 
manner. The as variation has the largest effect, ranging from 10% to 20% depending on the number of 
jets and Z?-jets. 

The WtZ background is important in the tetralepton channel. An uncertainty of 10% is assigned to 
the cross section, coming from the variation of renormalisation and factorisation scales. An additional 
uncertainty arises from the modelling of the additional jet in WtZ events. It is evaluated by varying 
parameters in simulation as described above. 


Other prompt lepton backgrounds: Uncertainties of 20% are assigned to the normalisations of the 
WH and ZH processes, based on calculations from Ref. [78]. An uncertainty of 100% is considered for 
triboson and same-sign WW processes. 
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Misidentified lepton charge background: This affects mainly the 2e-SS and efi-SS regions. Uncer¬ 
tainties on it arise from a variety of statistical and systematic effects. The main uncertainty comes from 
the limited statistical precision in the measurement of charge misidentification rates. These are treated 
as correlated among all px and I 77 I bins, but uncorrelated between the 2e-SS and efi-SS regions. This 
approach is comparable to treating each binned rate measurement as uncorrelated. Additional systematic 
uncertainties arise from the background subtraction in the Z ^ ee sample used to measure charge misid¬ 
entification rates and from the difference in charge misidentification rates between Z —> ee and ti events. 
The latter are evaluated by measuring charge misidentification rates in Z —> ee MC events, assigning them 
to ti simulated events and comparing the prediction to the number of true charge-misidentified events in 
the it sample. This uncertainty is found to be 10%. 


Fake lepton background: This is important in the same-sign dilepton, trilepton, and tetralepton chan¬ 
nels. The uncertainty on this background is estimated by propagating the statistical uncertainty on the 
measurement of the fake lepton efficiencies. Additionally, the normalisation of backgrounds in the con¬ 
trol regions (from prompt leptons or charge-misidentified electrons) is varied to estimate a systematic 
uncertainty on these efficiencies. The variation assigned (20-25%) depends on the composition of each 
control region, and is chosen to conservatively cover the largest uncertainty on the backgrounds. In the 
trilepton channel an additional uncertainty is considered by measuring the rates of real and fake leptons 
in two orthogonal regions, one with three or more jets and the other with one or two jets. In the same-sign 
dilepton channel, statistical uncertainties are dominant and no further systematic uncertainties are con¬ 
sidered. All uncertainties associated with fake leptons are considered to be uncorrelated among analysis 
channels and regions. 

In the opposite-sign dilepton channel the fake lepton background is small compared to the other back¬ 
ground contributions. An uncertainty of 50% is assigned to the fake lepton yield across all regions in this 
channel to cover the maximum difference between yields obtained from the simulation and from same- 
sign dilepton events in data. An additional uncertainty is assigned to cover the difference in shape of the 
distribution of the scalar sum of the transverse momenta of all reconstructed jets and charged leptons in 
the simulated and same-sign data events. 


7 Results 

The observed yields in the 15 signal and 5 control regions are shown together with the expected numbers 
of events in Table 4. 

The production cross sections atm and atu are determined simultaneously using a binned maximum- 
likelihood fit over all regions and discriminant bins considered in the analysis. The fit is based on the 
profile likelihood technique, in which the systematic uncertainties are treated as nuisance parameters 
with prior uncertainties that can be constrained by the fit. The calculation of confidence intervals and 
hypothesis testing is performed using a modified frequentist method as implemented in RooStats [79, 
80]. Significance is calculated using the asymptotic formula of Ref. [81]. 

A summary of the fit to all four channels with their corresponding fit regions used to measure the ffVT and 
tfZ production cross sections is shown in Figure 10. The normalisations of the ITZ and ZZ processes are 
determined as well, as described in Sections 5.3 and 5.4. 
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Figure 10; Expected yields after the fit compared to data in the five control regions (CR), used to constrain the tt, 
Z, ZZ and WZ backgrounds, the three signal regions in the opposite-sign dilepton channel (2LOS), the three signal 
regions in the same-sign dilepton channel (2LSS), the four signal regions in the trilepton channel (3L) and the five 
signal regions in the tetralepton channel (4L). In the two dilepton channels the fit also includes shape information. 
The “Rare SM” background summarises all other backgrounds described in Section 3 and mainly consists of the 
ttH, tZ and WtZ processes, which are the largest contributions to this background category in the dilepton, trilepton 
and tetralepton channels, respectively. The hatched area corresponds to the total uncertainty on the predicted yields. 


Table 7 provides a breakdown of the total uncertainties on the measured ttW and ttZ cross sections, 
determined by fitting each signal individually with the other fixed to its expected SM value. For both 
processes, the precision of the measurement is dominated by statistical uncertainties. For the ttW fit, 
the dominant systematic uncertainty source is the modelling of fake leptons and background processes 
with misidentified charge. For the ttZ fit, the dominant systematic uncertainty source is the modelling of 
backgrounds from simulation. 


Uncertainty 

O' ttW 

O'ttZ 

Luminosity 

3.2% 

4.6% 

Reconstructed objects 

3.7% 

7.4% 

Backgrounds from simulation 

5.8% 

8.0% 

Fake leptons and charge misID 

7.5% 

3.0% 

Signal modelling 

1.8% 

4.5% 

Total systematic 

12% 

13% 

Statistical 

+24% / -21% 

+30% / -27% 

Total 

+21% / -24% 

+33% / -29% 


Table 7: Breakdown of uncertainties on the measured cross sections of the ttW and ttZ processes from individual 
fits. Systematic uncertainties are symmetrised. 

The sensitivity to the ttW process is dominated by the same-sign dilepton channel, while the ttZ process 
is mainly measured in the trilepton and tetralepton channels. The result of the simultaneous fit of the ttW 
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Channel 

ttW significance 
Expected Observed 

tiZ significance 
Expected Observed 

leos 

0.4 

0.1 

1.4 

1.1 

2fSS 

2.8 

5.0 

- 

- 


1.4 

1.0 

3.7 

3.3 

M 

- 

- 

2.0 

2.4 

Combined 

3.2 

5.0 

4.5 

4.2 


Table 8; Expected and observed signal significances for the ttW and ttZ processes determined from the fit to the 
separate channels and from the combined fit to all channels. The significance for each signal process is calculated 
assuming the null hypothesis for the process in question and treating the other as a free parameter in the fit. 

and ttZ processes using all four channels is summarised in Table 8. The observed (expected) signific¬ 
ance of the measurements are 5.0cr (3.2cr) for the tfW process and 4.2cr (4.50") for the ttZ process. The 
background-only hypothesis with neither ttZ nor ttW production is excluded at V.lcr (5.90"). 

The result of the combined simultaneous fit to the two parameters of interest is 

cTttw = 369!^^ (stat.) ± 44 (syst.) fb = 369!^°° fb (4) 

and 

o'ttz - 176^43 ± 24 (syst.)fb = 176^g2 (5) 

Figure 11 provides a comparison of these measurements with NLO QCD theoretical calculations using 
M adgraph5_aMC @ NLO. 


8 Conclusion 

Measurements of the production cross sections of a top quark pair in association with a IT or Z boson 
using 20.3 fb“^ of data collected by the ATLAS detector in = 8 TeV pp collisions at the LHC have 
been presented. Final states with two, three or four charged leptons are analysed. From a simultaneous 
fit to 15 signal regions and 5 control regions, the ttW and ttZ production cross sections are measured 
to be (Ttiw - 369;^;gj°fb and cTttz - The fit to the data considering both signal processes 

simultaneously yields significances of 5.0cr and 4.2cr over the background-only hypothesis for the tfW and 
ttZ processes, respectively. All measurements are consistent with the NLO QCD theoretical calculations 
for ttW and ttZ processes. 
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uncertainty contours. The shaded areas correspond to 14% uncertainty, which includes renormalisation and factor¬ 
isation scale uncertainties as well as PDF uncertainties including as variations. 
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